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Abstract 

We study a neutron star with a quark matter core under extremely strong magnetic fields. We 
investigate the possibility of an Urea process as a mechanism for the cooling of such a star. We 
found that apart from very particular cases, the Urea process cannot occur. We also study the 
stability of zero sound modes under the same conditions. We derive limits for the coupling constant 
of an effective theory, in order the zero sound to be undamped. We show that zero sound modes 
can help kinematically to facilitate a cooling process. 
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I. INTRODUCTION 



The cooling of neutron stars is a fascinating subject that is far from being well understood. 
There are several mechanisms that contribute to the cooling of a star. Generally, neutron 
stars cool down by emission of neutrinos from the whole volume, or photons from the surface. 
In particular, neutrons stars with nuclear matter cool down initially due to neutrino emission 
via the so-called modified Urea process, where the emissivity scales as temperature T^, and 
later on, when the temperature becomes sufficiently small, due to photon emission from the 
surface of the star, that scales as T^. In the case of nuclear matter, the direct Urea process is 
dnematically forbidden, because energy and momentum cannot be conserved simultaneously 

ii,y,y,Q,y,y,0,y,3. 

On the other hand, neutron stars with quark matter, can facilitate the direct Urea process. 



This was demonstrated by Iwamoto 1^, The emissivity of the direct Urea scales as T^, 



which makes quark stars to cool faster than the corresponding ones with nuclear matter. 
The calculation of Iwamoto assumes that quark matter is ga ple ss. Direct Urea can also 
occur in dense nuclear matter with nonzero hyperon density 12 1, pion condensation jisj ]. 



and kaon condensation 



14| . Direct Urea can occur in all phases of quark matter, as long as 



the quarks are gapless. In sufficiently dense quark matter, quarks might pair in a nontrivial 



way, forming Cooper pairs. This is the phenomenon of Color Superconductivity [15|, |16 |. 
In this case, fermionic excitations might acquire an energy gap, thus making impossible the 
direct Urea process for phases like the Color-Flavor-Locking (CFL) 17|]. However in several 
cases, some fermionic excitations might remain gapless, and if kinematically allowed, they 
can facilitate a direct Urea process 18 1. 

The case of neutron stars with strong magnetic fields is an interesting one. Although 
weak magnetic fields do not seem to infiuence a lot the cooling process, strong magnetic 
fields might alter the picture completely. This is because charged particles populate Landau 
levels. Under strong magnetic fields, where eB ~ /i^ (/i being the characteristic chemical 
potential), only the very first Landau levels are occupied. So far, several neutron stars have 
been observed with magnetic fields of 10^^-10^^ Gauss. It is possible that the magnetic fields 
inside the star could be even higher. The emissivity of the direct Urea process for the case 



of a neutron star made of nuclear matter with strong magnetic fields was studied in 19|, |20 |. 
It was shown that for eB > /x^, direct Urea reactions are open for an arbitrary proton 



9. 



concentration. Due to the fact that direct Urea has in general a larger emissivity compared 
to modified Urea, neutron stars with strong magnetic fields will cool much faster than stars 
without. By comparison of observational data, constraints can be imposed regarding the 
strength of the magnetic fields in the neutron stars, or vice versa, predictions can be made 
for anomalously cold neutron stars. 

Although the case of neutron stars made of nuclear matter with strong magnetic fields 
has been studied and is somewhat understood, the case where the neutron star is made 
(partially) of quark matter under strong magnetic fields has not been studied to the same 
extend. Color superconducting quark matter under high magnetic fields has been studied 



in 211]. If quark matter exists inside neutron stars, in order to fill only the lowest Landau 
level, magnetic fields qB ~ /i^ are required. For a quark chemical potential of the order of 
100 MeV, which is expected for quark matter in a neutron star, the corresponding magnetic 
field is of the order of 10^^ Gauss. As we already mentioned, the highest magnetic fields so far 
observed in magnetars are 10^^ — 10^^ Gauss. However, it is probable that stronger magnetic 
fields develop in their interior. From this point of view, the study of strong magnetic fields 
inside a neutron star with quark matter becomes an interesting problem. 

In this paper we shall argue that the direct Urea does not occur in general in this case. 
In fact we shall show that at extremely high magnetic fields, the "triangle inequalities" , i.e. 
the simultaneous conservation of energy and momentum is not possible and therefore direct 
Urea is not open, with the exception of very specific choices of the chemical potentials. 

This paper is organized as follows: In section II we show why direct Urea is not open 
in general for quark matter in strong magnetic fields. In addition, we point out under 
what specific conditions such a reaction can take place. In section III we study the zero 
sound modes of quark matter under strong magnetic fields and we show how they can help 
in satisfying the triangle inequalities, making possible a cooling reaction. We conclude in 
section IV. 

II. THE IMPOSSIBILITY OF DIRECT URCA 

The direct Urea reaction in the case of quark matter consists of two reactions 

M + e — > c? + z/, (1) 



d ^ u + e + u. 



(2) 



A similar reaction with a strange quark instead of a down, can also take place. In the 
presence of a magnetic field, charged particles fill the Landau levels. The wavefunction of a 
charged particle is 

= /r T \V2 exp[-«^t.t + ipyV + ipzz]u 



(3) 



where 



I {E^ + m)I^,p^{x) \ 


Pzlu,py{x) 



and 



u 



1 



[2E,(E, + m)]V2 





{E^ + m)J, 



\ 



X] 



[X) 



The symbols | and | denote up and down spin states. We have chosen the magnetic field B 



to be along the z axis. The symbol v indicates the Landau level, running from 0,1,2, 



and Lz are nonphysical normalization lengths for the y and z axis, and they should drop 
out at the end of the day. The function ly^p^ is given by 



1/4 



1 



qB\ 

TT ) (2^Z/!)l/2 



exp 



IqB (x-^ 



{.BY'Hx-^ 



(4) 



where Hy is the well known Hermite polynomial. We have assumed that the electric charge 
q is positive, but similar formulas can be written down for a negative electric charge. The 
energy eigenvalues are Ey = (p^ + + 2z/gi?)^/^, where m is the mass of the particle. Here 
we are going to assume that eB ~ /i^ and therefore all the charged particles occupy the 
lowest Landau level, i.e. u = 0. If the quark matter is in /? equilibrium. 



(5) 



In order for ([T]) to take place, momentum along the z axis and energy should be conserved. 
This leads to 



Pid + fnl + Ey for energy. 



(6) 



Pzu + Pze = Pzd + Pi^ for z momentum. (7) 

Similar equations hold for the case of the strange quark, where we can substitute c? — s> s in 
Eqs. ([H [21 El ED • Obviously, the /5 reactions with the strange quark are somewhat suppressed 
due to the Cabibbo angle. Since we have assumed that all the charged particles occupy only 
the lowest Landau level, this means that the particles fill all the states between —pp to 



Pf, where Pf = y — Having chosen the magnetic field along the z axis, up quarks 
in the lowest Landau level have their spin aligned along the positive z axis, whereas down 
and strange quarks, as well as electrons have their corresponding spins along the —z. The 
electroweak sector involves only left handed particles, and therefore at the limit where the 
mass of the particles goes to zero, only the up quarks with momentum —pFu, down quarks 
with momentum pFd, and electrons with pFe would participate in /5 reactions. However, 
since the particles do have a mass, up quarks with +Pfu, as well as down quarks with 
—pFd and electrons with —pFe can participate, although the amplitude will be suppressed as 
m/fi, with m and /i being the corresponding masses and chemical potentials of the particles. 
Similarly for the reaction ([21) we have the following equations from energy and momentum 
conservation 



VpL + ml + Vp1e + ml + Ei, = ^pI^ + ml for energy, (8) 

Pzu + Pze +Pp = Pzd for z momentum. (9) 

Now we shall argue that the set of Eqs. ([SI H [HI IH]) cannot be satisfied in general. Firstly, 
let's consider the case where up, down and electrons have velocities aligned along the +z 
axis. Quark matter if present inside neutron stars, is expected to be at high densities, which 
means that m « fi and m << pf- We can expand Eqs. ([6l [HD around the corresponding 
Fermi momenta and we get respectively 

2 2 2 

m,, mi mi „ , 

PFu+7r^+PFe + 77^=PFd+7r^+ E„ 10 
2/X„ 2/ie 2/id 

2 2 2 

mf, mt mi , , 

PFn + TT^ + PFe + 77^ + = PFd + 77^- (H 
2/i„ 2/ie 2/irf 

Subtracting ([7]) from ( ITOl) and ([9]) from ( ITTl) = pf^ , we get 

E,-p, = -A, (12) 

and 

Eu-Pv = A, (13) 



where 



^_ ml ml ml 



2fid 2flu 2/Xe 

We can easily see now that the Urea process cannot occur in this case. First of all, since 
the energy of neutrinos is ~ T, the above equations require also A ~ T, which is not true in 
general. However even if this is true, since Ei,^p > pj,^p, Eqs. ( IT2l) . and ( |T3l) cannot be satisfied 
simultaneously provided A ^ 0. Now we can consider a case where one of the particles has 
velocity with opposite sign with respect to the rest. The most important case is to consider 
the up quark with p^u = —ppu- In this case there is no m//i suppression of the amplitude. 
Because of the flip of the sign of the momentum, Eqs. f|T2l [T3|) become 



E^-p^ + A = 2pFu, (14) 

Ep-pp + 2pFu = A. (15) 

This implies that the reactions ([H [2]) can take place if 2ppu = ^- This relation can be 
rewritten if we expand around /x„ as 

2 2 2 

mt mt ,^ . 

2/i„ = ^ + ^ ^. 16) 

2/irf 2/i„ 2/ie 

Although in principle such a relation is not forbidden, it is highly unlike to be satisfied, 
as it requires very low chemical potential for the up quark. This can be improved by 
considering strange quarks instead of down, because the mass of the strange quark is much 
larger compared to the down. In that case, 

where we omitted the subleading terms of ( |T6l) . If it is the electron and not the up quark 
that flips velocity, we can carry the same analysis, concluding /le = m^/(4/is). However, in 
this case the reaction will be suppressed by factors of {rriumf,)'^ / (fiuf^e)"^ ■ Eq- (pTj) (or the 
corresponding with e instead of u) is very restrictive and hard to be satisfied. It means that 
only for a very specific value of the chemical potential, the Urea process can occur. 

We can now impose overall electric neutrality for the quark matter. This implies that 
{2/3)nu = {l/3)nd + {l/3)ns + rig, where n stands for the number density of the particular 
particle. To calculate the number density we have to make the substitution 

dpxdpy qB 
47r2 ^ 2tt' 



Upon making this substitution, the number density is 



qB , . 

n = —^PF. (18) 



After expanding for small masses, the neutrality condition becomes 

mi rn?, rrii 2rn?, . . 

/i, + /irf + 9/i, - 4/i„ = + + 9-^ ^. 19 

2/ie /i„ 

If we assume that the reactions of ([1]), and ([2]) take place for the strange quark, thermal 
equilibrium suggests /i^ = fiu + f^e- As we showed, the same reactions cannot hold simul- 
taneously for the down quarks, because it would imply a chemical potential for the quarks 
of the order of or less. Therefore f3 equilibrium between up and down quarks cannot 
be established through Eqs. ([HE]). In this case, we have four parameters /i^, /i^, and 
/ie, where one is fixed by Eq. (fT7|) . one by chemical equilibrium between strange, up and 
electrons, and one by electric neutrality. There is one left free parameter, and in order for 
direct Urea to occur, the neutrality condition (after using the chemical equilibrium condi- 
tion, Eq. ( |T7|) . and having kept only the first term from the right hand side of Eq. ( fT9l) ). 
reads 10/is + fid = 15mg/(4/is). If for any reason the down quarks are in equilibrium with 



the other species, fig = fid- In this case, the above condition becomes /i^ = msA/l5/44. 
This is of no interest, because direct Urea occurs only for a single value of /i^, which is 
in addition unrealistically low. Furthermore, in this case the expansion we have done is 
not valid any more since /i^ is smaller than m^. If the triangle inequalities are satisfied for 
the electron instead for the up quark, i.e. = m^/(4yUs), the condition for direct Urea is 



Hd = 'ijJ's — ll"^s/(4Ai«)- Again, in case where /i^ = /i^, the condition reads /is = y/ll/S. 
As we mentioned, it is not expected that /i^ = /i^, since chemical equilibrium between up, 
down and electrons cannot be established in general. Our conclusion is that the direct Urea 
process can occur only if Eq. (fTTl) is satisfied (with either the up quark or the electron) and 
fid is such that satisfies the neutrality condition. 



III. ZERO SOUND MODES 



As we argued, in the case of quark matter at extremely strong magnetic fields, the direct 
Urea process cannot occur in general. Only under very specific conditions and probably 
unrealistic, such a process can take place. In the case of nuclear matter (without magnetic 



7 



fields), the direct Urea does not occur either (unless there is nonzero hyperon density or 
pion, or kaon condensation), because the triangle inequalities are not satisfied. Instead, the 
modified Urea, with the presence of a bystander particle takes place. In our consideration, a 
bystander particle cannot help. This is because even the extra particle, either it is up, down, 
strange or electron, it should also be aligned along the z-axis, and therefore it is difficult 
again to satisfy the triangle inequalities. As in the direct Urea, perhaps such a process can 
occur, but again only under very specific values of the chemical potentials. 

In this section we are going to study the zero sound in this particular setup, and we 
are going to investigate if such modes can help establish a cooling process by ensuring 
that the triangle inequalities are satisfied. First let's consider how the zero sound modes, if 
undamped, can facilitate simultaneous conservation of energy and momentum. Let's assume 
for the moment that the zero sound modes are undamped and that their velocity is fo < 1. 
We can consider the following reactions (which are modifications of the (3 reactions ([H [2]) 
with the inclusion of a zero sound) 

u + e ^ d + z/ + zs, (20) 

(i + zs— >M + e + j^, (21) 

where zs stands for zero sound. The equations that correspond to the previous Eqs. f[T^ [T^ 
are respectively 

E,-p,-{l-vo)k = -A, (22) 

and 

E^-p^ = A-{l-VQ)k. (23) 

k is the momentum of the zero sound, and the term —{l—VQ)k is uo — k = v^k — k = —{l—Vo)k, 
where we used the dispersion u = v^k for the zero sound. From Eqs. (l22l 1231) . we see that 
energy and momentum can be conserved, if we require (1 — VQ)k = A. In some sense, 
the above mechanism corresponds to a "modified Urea" with the bystander being the zero 
sound, rather than a particle. We have already mentioned that modified Urea with a charged 
particle as a bystander would not help, as it would face the same limited phase space with 
direct Urea due to the fact that the bystander will have to be aligned along the z axis 
with a Fermi momentum. Here, we see that in principle, zero sound modes can make this 
"modified" version of the Urea process to occur at high magnetic fields. Of course for such a 
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process to take place, the zero sound modes should be undamped. In the rest of this section 
we are going to study the zero sound, in the case of quark matter under strong magnetic 
fields, using an effective Lagrangian. 

The zero sound is a collective mode in Fermi liquids, involving particles near the Fermi 
surface. It was predicted by Landau 2^ in 1959 and it was observed experimentally in liquid 



He^. The relativistic extension of the Landau Fermi liquid theory (applicable for example 



at high density QCD) was developed by Baym and Chin 23|]. In that paper, the zero sound 
velocity was calculated for a relativistic Fermi liquid using different effective models. Here, 
we calculate the zero sound velocity for a relativistic Fermi liquid that is under ultra strong 
magnetic fields, and therefore all charged particles fill only the lowest Landau level. 



The Boltzmann equation that governs the number density in the coUisionless case is 24 1 



dn de dn de dn ^ . ^. 

dt dp dr dr dp ' 

where e and n are the energy and density distribution of the quasiparticle. We consider 
small deviations from the equilibrium distribution 

fi(p,r,t) =no(p) + 5n(p,r,t). (25) 

The infinitesimal change in the energy is 

Se = ti' [ fSh'dr', (26) 



where J dr' is the integral over the phase space, the trace is over spin states, and / is the 
usual Landau interaction function. If we expand Eq. ( 124|) and keep only up to first order 
variations, we have 

d5n ^ deo d5n _ d5e dno _ ^ 
dt dp dr dr dp ' 

with eo and uq being the energy and number density equilibrium distributions. Because all 

the charged particles are in the lowest Landau level, the equilibrium distributions for the 

density and the energy depend only on the z component, deo /dp = v, where v = vpn, n 

being the unit vector in the direction of p. It is understood that n can be along the ±z 

axis. On the other hand, dn^/dp = —v6{e — ep). We are seeking wave solutions of the form 

Sfi = 6{e - eF)0{n)e'^^-''-'^'\ (28) 



Before we investigate the solution of Eq. (!27|) . we derive the Fermi interaction function 
/, as it is quite different from the case where the magnetic field is absent. Let's assume that 
the Lagrangian includes an interaction term of the form 

I^int = ig^Y^A^, (29) 

where A is a massive vector field. In this case we can write the amplitude of two fermions 
scattering by exchanging A in the space coordinates as 

Tfi = ig"^ j d^xd'^x'ffi{x)D{x - x')j^fi{x'), (30) 

where jj^ = ijj'-f^ip and the vector propagator D is 

r d\ expHg.(a:-xO] 
^^"^ 7(2^ ' ^^^^ 

where M is the mass of the vector boson. We can focus on the case of a positive fermion 

(for example the up quark), where the wavefunction is given by Eq. (3), with v = Q for 

the spinor, since we consider only the lowest Landau level. Let's denote by the subscripts 1 

and 2 the incoming particles and 3, and 4 the outgoing ones. Eq. ( l30l) involves integration 

over d'^xd^x'd'^q. Because the spinor of the wavefunction in the lowest Landau level depends 

only on x, we can easily perform the integration over dt'dy'dz' that gives the product of 

delta functions {2'7iYLyLz6{qo — E2 + E4)5{qy — Py2 + Pyi)S{qz — Pz2 + Vza)- Integration over 

dq^dqydqz is now trivial due to the delta functions, and finally after integrating over dtdydz, 

we get 

dxdx'^^-^B,{x)B'^{x'), (32) 

where B^{x) = u\x)'-f'^u\x), and the delta functions denote overall energy and mo- 
mentum {y and z axis) conservation. The Landau interaction function is defined as 
Tfi = i{27TY6{E)6{py)6{pz) f ■ There are two distinct channels, we should consider. One 
is the direct forward scattering where Ei = E^ and E2 = E4 (with similar relations for the 
corresponding momenta), and the second is when Ei = E4 and E2 = E3 that corresponds to 
an energy and momentum exchange between the two particles. For the first case the integral 

f dq^ 1 exp[-M|x-x'|] . . 

J 27^^3^^"Pt^^^^" - " = 2M ■ 

The Landau interaction function is 

(34) 

in 



For convenience prime energy and momentum correspond to the ones with subscript 2 (in 
the direct amphtude subscript 4 as well). We proceed now to calculate d6e/dr as we need 
it for Eq. ( 1271) . The spin trace of Eq. ( l26l) is trivial as all particles have spin up. Plugging 
in the solution of Eq. fl25]l , we have 

or J ZTT ZTT 

After integrating over dpy'dxdx' we get 

f S^/<'-^(-'--.)KnO(l-||;). (36) 

There are only two directions for n', either up or down. After we perform the integration, 

We can now rewrite Eq. (1271) using the solution (1281) as 

qB ( ^ „ ^ ^ 2E'^p — 



This equation leads to a homogeneous system of two equations 

[uj - VFk{l + x)]z>(+) - Vfkx5v{-) = 0, (39) 
[uj + vpkil + x)]z/(-) + ?;/A;x(5i>(+) = 0, (40) 

where 

qB 



X 



47r2 M^EfPf 

S = {2E'^—m?)/m?, and ± denotes the direction of the momentum p along or opposite to the 
direction of k. The latter has two distinct directions along the ± z axis. The determinant 
of this system of equations should be zero if a nontrivial solution exists. This gives the 
following dispersion relation 

cu2 = 4A;2[1 + (1-5V + 2x]. (41) 

Generally, in order for the zero sound to be undamped, its velocity should satisfy vq > vp. 
This condition can be easily proven upon making the requirement that the frequency uo has 
no imaginary part. However, in the case of a strong magnetic field, this condition is different. 
For u) to be real, Eq. ( 14T1) gives 

X < (42) 
11 



Assuming m << Ep, the above equation becomes 

q M 

— < 43) 

Unless the couphng constant satisfies the above equation, the zero sound modes will be 
damped. We should note that if we were to impose the general condition vq > vp, the upper 
bound for the coupling g would have been 

g Mm 



2n ^/2^Ef 

which is suppressed compared to the one we derived above, by an extra factor of m/EF- 
There are two points that we would like to stress here. The first one is the difference between 
our constraint about the undamped zero sound and the Landau general one {vq > vp) that 
holds in the case without magnetic field. The Landau constraint from the mathematical 
point of view comes from the requirement that uj has no imaginary part, which refiects the 
simple physical fact that zero sound modes with velocity less than vp lie in the continuum, 
and they are unstable as they can decay back to a particle and a hole. With vq > vp, such a 
situation is avoided as the zero sound mode cannot decay to a particle-hole pair preserving 
energy and momentum simultaneously. The case we study here is qualitatively different. 
Because the particles fill only the lowest Landau level, the problem is one dimensional 
rather than three. The argument we used above, cannot be used now: there can be zero 
sound modes with velocity smaller than vp that cannot decay to particle- hole pair, because 
there is no continuum in this case. In the three dimensional case u = vpk cos 6, with 6 being 
the angle between the vf and k vectors. For a proper choice of 9, u can take all the values 
between and = |vF||k|. In our case though, because the problem is one dimensional, 
fermionic excitations above the Fermi surface can lie on the line uo = vpk, but not below 
(or above). Therefore zero sound modes with vq < vp, cannot decay to a particle-hole pair. 
Our constraint of Eq. (42) or (43), show that there is a dynamical instability if it is not 
fulfilled, but it is not related to Landau damping. 

The second point has to do with the notion of zero sound in one-dimensional systems. 
Generally, zero sound modes are viewed as distortions of the Fermi surface along the forward 
scattering direction, while first sound modes correspond to simple radial oscillations of the 
Fermi surface. It is obvious that such a distinction is not valid in the case we investigate, 
as the Fermi sphere is transformed to an one-dimensional column. From this point of view, 

19 



the distinction between zero and first sound modes becomes fuzzy. However, we can still 
talk of zero or first sound with respect to the quantity ut, that is the energy of the mode 
times the mean time between collisions for the quasi-particle. ur >> 1, meaning that the 
frequency of the mode is much larger than the frequency of collisions, corresponds to zero 
sound, whereas ujt << 1 correspond to the case where the hydrodynamic limit is valid and 
the modes correspond to first sound. Generally, zero sound exists in Fermi liquids at low 
temperatures since r ~ T"^. However in our case, because of the difference in the phase 
space, we should examine how r scales. This is easy to do if we recall that the particles 
that can interact with each other are the ones that have momenta within T from the Fermi 
surface. This means that the effective density of interacting particles is An = qBT/n'^. The 
effective cross section for a collision is a ~ {qB)~^, and the mean free path (and consequently 
r) should scale as / ~ T~^. For energy u larger than a few keV (which is a typical neutron 
star temperature), it is safe to say that ujt » 1 and therefore we indeed have a zero sound 
mode. 

We should also stress here, that the quasi-particle energy is given by 

off, \ / \ dpy' dpz' , , ,s 

ep = ep+ / fiPy^Pz;Py'^Pz')n{Pz')- —■ (44) 

J Zn Zn 

Plugging the expression of / for the direct amplitude and keeping in mind that to lowest 
order —pp < Pz' < Pf, we get 



^P = 4 + T72^^ (45) 



1^ 

where n is defined in Eq. (1181) . This ag rees with the corresponding relativistic Fermi liquid 
(and zero magnetic field) formula of 23|] and with the result of 25|. We should also mention 
that the zero sound velocity vo should be obviously smaller than the speed of light. It turns 
out that if 5^(1 — Vp) > 1, vq is smaller than one for any value of x or correspondingly of 
the coupling constant g. If we choose a typical chemical potential of 500 MeV for a quark 
of mass 10 MeV (say down for example) 5^(1 — Vp) ~ 10^ > 1 which means that t;o < 1 no 
matter what the coupling constant g is. 

Now we discuss the exchange amplitude of the two interacting particles. In this case 
El = and E2 = E^. Similar relations hold also for the momenta along the z and y axis. 
The integral that corresponds to Eq. ( |33l) now becomes, 

dqx I , . I exp[-C|x-x'|] 
n 7TT7 exp — iOj, X — X = : , 46 
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FIG. 1: Left panel: 5 as a function of the mass M (in MeV) for three different chemical potentials 
and magnetic fields, namely fi = 400 MeV (thick solid line), /i = 500 MeV (thin solid line), and 
H = 600 MeV (dashed line). In all three cases, the magnetic fields have been set to \J2qB = /i, 
which is the lower possible magnetic field in order to populate only the lowest Landau level. Right 
panel: 5 as a function of M for a fixed /i = 500 MeV, for three different magnetic fields, namely 
^/qB = 500/^2 MeV (thin solid line), 500 MeV (dashed line), and 600 MeV (thick sohd line). 



where C = + Q'^ + ~ QxtY^"^- The Landau interacting function is 



f{Pv,Pz;Py',Pz 



g^qB 
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dxdx'- 



-({x-x'l 



m 



c 



where 



(47) 



K{x,x\py,py,) = -^[(x - + {x- + (x' - + {x' - (48) 



As in the previous case, we are looking for solutions of the form of Eq. fl28|l for the Boltzmann 
Eq. (|27l) . The dispersion relation reads 



(a; — v.k)i>(n) = v.k 



f dpy \^f(^p^^p^.p^,^ppn)u{n) + f{py,Pz]Py', -pFn)u{-n)] ) . 



Pf \J 47r2 



(49) 

By using Eqs. (147 1 l48l l49l) . we can perform the integrations over dpy/dxdx'. Eq. (l49l) reads 



(cj =F VFk)i){n) = ±VFk{fii){n) + /2Z>(-ra)), 



(50) 



where 



ha 



m 



4v^7r3/2 EfPf 



,C.lJ{2qB) 



C: 



-Erfc 



1,2 



Cl,2 



(51) 
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FIG. 2: Left panel: An upper bound on g/{2Tr) as a function of the mass M (in order vq < 1) for 
the parameters chosen in the left panel of Fig. 1. Right panel: An upper bound on g/{2TT) as a 
function of the mass M (in order vq < 1) for the chosen parameters of the right panel of Fig. 1. 

As we have mentioned ( = {Qy + + M"^ — q^Y^"^- Q'o = (since both the colliding particles 
are on the Fermi surface). It is easy to show that the integral of Eq. fl47p does not depend on 
Py, and therefore we can set it to zero. On the other hand, the integral over dpyi is dominated 
by the values of pyi close to Py. For this reason, and in order to obtain an analytic result, we 
have neglected the Pyi dependence of C. In the case where the colliding particles are moving 
in opposite directions, = '^Pf » %■ Therefore this approximation is justified. If the 
colliding particles are moving towards the same direction, = 0, and this approximation 
is justified only if << M. According to this, (i = M (for same direction collision) and 
C2 = {4^p'f + M'^Y^'^ (for opposite direction collision). The dispersion relation is given again 
by Eq. (HTj) . where in this case x = fi and 6 = /2//i- However unlike the previous case 
where 6 » 1, here 5 < 1 as it can be seen in Fig. 1. Due to this, as long as 5 < 1, the zero 
sound waves (provided there is no direct amplitude) are undamped. Again the zero sound 
velocity should he vq < 1. This condition leads to a constraint in g that is depicted in Fig. 2. 

In principle, instead of deriving the dispersion relations separately for the direct and 
exchange amplitudes, we should have added the two contributions / = fjjir + /exc- However, 
we investigated separately the two cases for two reasons. Firstly, there are interactions like 
QCD that do not have contributions from the direct amplitude due to the overall gauge 
structure. For example, in QCD the direct amplitude is zero because it will be proportional 
to the trace of a single Gell-Mann matrix which is zero. On the other hand, in case where 
both amplitudes are nonzero, by inspection of the x below Eq. ( I40i) (let's call it xoir) and 
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FIG. 3: The ratio xoir/ fi as a function of M for the set of parameters chosen in the left panel of 
Fig. 1. 

/i, we can conclude that for most of the cases xdit > fi. This is depicted in Fig. 3. In fact 
we can see that for low M the ratio becomes large and therefore the exchange amplitude 
is small compared to the direct one. In this case Eqs. (42, 43) still hold. For larger M, 
fi ~ XDir- However, in any case Soir >> ^Exc- Therefore Eq. (42) should be rewritten in 
the general case as 



where a = 1 + fi/xoir- In Fig. 3 we see that 1 < a < 2. As before, for 6 » a, or equivalently 
for m << Ep, the constraint reduces back to Eq. (43). It is safe to say that when both 
present, the direct amplitude dominates over the exchange one. 

As it can be seen in Fig. 3, for large M, XDir/ fi — ^ 1- This is expected, because large M 
means small range for the interaction, and therefore direct and exchange amplitudes become 
similar when momenta are aligned, since the interaction tends to be very localized (due to 
large M). The same is not true if the interacting particles have opposite momenta. One 
can easily check that direct and exchange amplitudes behave differently in this case even at 
large M, because of the different Dirac structure. At the chiral limit m = 0, the exchange 
amplitude is zero, because the vector boson cannot flip the chirality, whereas in the direct 
amplitude there is no need for flipping the chirality. In fact at the chiral limit all amplitudes 

Ifi 



except the direct one with opposite momenta for the interacting particles are zero. The 
reason why for example the direct amplitude with aligned momenta Pz = p'^ is zero at the 
chiral limit whereas for pz = — is not, has to do with the simple fact that particles move 
with the speed of light and due to the localized interacting force, they can never interact if 
Pz = p'z unless they are on top of each other. If pz = —p'z, they will always meet no matter 
how far they start from each other, provided their transverse distance is within 1/M. We 
should also mention that 6 in Fig. 1 saturates to 1 for larger M than the ones depicted in 
the figure. 

IV. DISCUSSION 

In this paper we studied if the direct Urea process can occur inside a neutron star with 
dense quark matter under extremely high magnetic fields. When qB ~ /i^, only the lowest 
Landau level is populated. As we showed, the lack of transverse momentum makes impossible 
to satisfy simultaneously energy and momentum, and apart from very specific cases that we 
pointed out, direct Urea cannot take place in general. In search for finding ways to save the 
Urea process, we investigated sound modes and in particular zero sound. We showed that 
in principle if such modes are undamped, they can help facilitating a variant of the Urea 
process, that in some sense is a type of Modified Urea. Furthermore, by using an effective 
model where fermions can interact via the exchange of a vector boson, we found under 
what constraints, we can have undamped zero sound. We derived the dispersion relations 
and we calculated the zero sound velocity. We found that unlike the "usual" 3-dimensional 
case (with absent magnetic field) where Vq > vp in order to have undamped waves, in the 
case studied here, vq can be even smaller than vp- This is due to the fact that there is no 
continuum in this case, since the Fermi sphere has collapsed to an 1-dimensional column, 
and therefore modes with velocity smaller than vp cannot dissolve back to a particle and a 
hole. 
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